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ABSTRACT 

A new approach for the classification of the cosmic web is presented. In extension of the 



previous work of Hahn et al. ( 2007| ) and |Forero-Romero et al. (2009! the new algorithm is 
based on the analysis of the velocity shear tensor rather than the gravitational tidal tensor. 
The procedure consists of the construction of the the shear tensor at each (grid) point in space 
and the evaluation of its three eigenvectors. A given point is classified to be either a void, 
sheet, filament or a knot according to the number of eigenvalues above a certain threshold, 0, 
1 , 2, or 3 respectively. The threshold is treated as a free parameter that defines the web. The 
algorithm has been applied to a dark matter only simulation of a box of side-length 64ft," 1 Mpc 
and N = 1024 3 particles with the framework of the WMAP5/ACDM model. The resulting 
velocity based cosmic web resolves structures down to < 0.1/i _1 Mpc scales, as opposed to 
the l/i~ 1 Mpc scale of the tidal based web. The under-dense regions are made of extended 
voids bisected by planar sheets, whose density is also below the mean. The over-dense regions 
are vastly dominated by the linear filaments and knots. The resolution achieved by the velocity 
based cosmic web provides a platform for studying the formation of halos and galaxies within 
the framework of the cosmic web. 



1 INTRODUCTION 



The visual appearance of the large scale structure (LSS) of the uni- 
verse is of a web, the so called cosmic web ( |Bond et al.|[T99"6) . 
The large scale distribution of galaxies in the observed universe as 
well as the distribution of the dark matter (DM) as inferred from 
its gravitational lensing and reconstructions from large galaxy sur- 
veys give the appearance of mass and light distributed in a web-like 
structure dominated by linear filaments and concentrated compact 
knots, thereby leaving behind vast extended regions of no or a few 



galaxies and of low density (Kitaura et al. 2009 


Jasche et al. 2010 


Munoz-Cuartas et al.||201lf Wang et al. 201 


i. Direct mapping 



of the mass distribution by weak lensing reveals a time evolving 
loose network of filaments, which connects rich clusters of galax- 
ies ( jMassey et al.|200"7) . The extreme low resolution of the weak 
lensing maps cannot reveal the full intricacy of the cosmic web, and 
in particular the difference between filaments and sheets, yet they 
reveal a web structure that serves as a gravitational scaffold into 
which gas can accumulate, and stars can be built. This poses an 
intriguing challenge of the mathematical classification and quan- 
tification of the cosmic web. The motivation for such an endeavor 
is twofold. On the one hand the cosmic web is there and so we want 
to describe it mathematically. On the other hand the web classifi- 
cation might provide some extra parameters that quantify the envi- 
ronment within which DM halos and galaxies form. The properties 
of galaxies are observed to depend on their environment (Dressier 



[TM)l|Blan"ton et al.|2005||Avila-Reese et al.|2005||Gao et al.|2005l 
Maulbetsch et al. 2007 , Forero-Romero et al. 201 1 ). Further more 



the dynamics of sub-halos and satellite galaxies suggests a possi- 
ble dependence on the environment within which their parent halos 
reside ( |Knebe et al.|2004||Libeskind et al.|2005|[20TTT >. A web clas- 
sification might extend our theoretical tools for understanding such 
a dependance. 

Translating the visual impression of a cosmic web into a math- 
ematical formulation is not trivial. The thorough review of |Colberg| 
|et al.|p008} of void finders provides some perspective to the gen- 
eral issue of web classifiers. This has been pursued along two dif- 
ferent lines, the geometric and the dynamic approaches. The geo- 
metric strategy focuses on the point process exhibited by the distri- 
bution of galaxies or DM halos in simulations, say, and describes it 
mathematically. This has been often applied independently of any 
dynamical context (Lemson & Kauffmann 199 9) |Novikov et al.| 
12006) |Arag6n-Calvo et al.||2007[ jSousbie et al.||2008|>. The dy- 
namic approach has its roots in the seminal work of |Zeldovich| 
{1970}, which led to the 'Russian school of structure formation' 
(e.g. |Amold et al] {[982); |Klypin & Shandarin| {1983) ) The Zel- 
dovich approximation has been the first analytical tool that enables 
the tracking of the formation of aspherical objects, thereby account- 
ing for the formation of voids, sheets, filament and knots. A recent 
application of that approach to the cosmic web is presented in |Lee| 
|&Lee|pQ08) 

The geometric approach uses objects (galaxies or DM halos) 
which are per definition at a position of locally high density. |Hahn| 
|et al.| {2007) proposed a very attractive dynamic classification of 
the cosmic web. They presented an algorithm by which each point 
in space can be classified as being either a void, sheet, filament 
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or a knot point. Namely, an attribute is assigned at each point in 
space much in the same way as it is characterized by its density 
and velocity. This attribute depends, in turn, on the the dynamics of 
the LSS, namely on the density and/or the velocity fields. 

The |Hahn et al.||2007} dynamic web classification is based 
on counting the number of positive eigenvalues of the tidal tensor, 
i.e. the Hessian of the gravitational potential. The number ranges 
over 0, 1, 2 and 3, corresponding to a void, sheet, filament and a 
knot attribute. This web classifier is of a dynamic nature and it is 
closely related to the equations of motion that dictate the dynamics 
of the growth of structure. One can consider the web classifier to be 
a dynamical field that characterizes the large scale structure of the 
universe, a field that accepts only four possible numbers. This ap- 
proach is driven by the Zeldovich approximation and its web classi- 
fication agrees remarkably well with the visual impression one has 
in viewing the large scale structure (LSS) on linear and quasi-linear 
scales, i.e. scales larger than a few Megaparsecs. However, in go- 
ing down to smaller scales the Hahn et al approach fails to recover 
the fine web exhibited by high resolution numerical simulations. A 
partial remedy to this shortcoming was provided by Forero-Romero 
|et al.| |2009 1 (FR09), who relaxed the assumption of the null thresh- 
old which is used to classify the web. FR09 argues that the value 
of the properly normalized threshold should be around unity. In- 
deed the cosmic web defined by a non-zero threshold provides a 
much better description on the non-linear cosmic web, down to the 
Megaparsec scale. 

Inspection of high resolution N-body simulations reveals a 
web-like structure on very small scales as well. In simulations 
aimed at resolving galactic scale halos the web is observed down 
to the virial radius of such halos (Libeskind et al. 2011). Study- 
ing such simulations we have tried to extend the approach of 
FR09 down to sub-Megaparsec scales, but have failed. Namely, the 
method has not been able to resolve the fine web of these scales. 

In the linear regime of the gravitational instability the velocity 
and gravitational fields are essentially identical, up to some scal- 
ing that depends on the cosmological parameters. Hence, the web 
classification of |Hahn et aklpOOT} and FR09 can be reformulated 
in terms of the (velocity) shear tensor. It follows that in the linear 
regime the tidal tensor based cosmic web and the velocity tensor 
based web are identical (hereafter these are defined as T-web and 
V-web, respectively). Going to the fully non-linear regime we ex- 
pect the T- and V-web to depart. Here, we study the cosmic web 
as revealed by studying the analysis of the fully non-linear velocity 
field. The algorithm proposed here is tested against a high resolu- 
tion DM-only N-body simulation. 

A general remark is due here, before proceeding to the main 
body of the paper. By itself the cosmic web is an ill-defined en- 
tity. The transition from knots to filaments, sheets and voids seems 
to be smooth and gradual, yet visual impression strongly suggests 
that the LSS is correctly characterized by a cosmic web. This res- 
onates somewhat with the issue of DM halos. Again visual impres- 
sion strongly suggests that the dynamics of gravitational instability 
leads to the formation of bound compact entities called DM halos, 
yet these objects have no strict boundaries and therefore their defi- 
nition would rely on somewhat arbitrary chosen free parameter(s). 
However, even if the gravitational collapse in an expanding uni- 
verse is still not fully understood a lot of insight and understand- 
ing is provided by the spherical top-hat collapse model (|Gunn &| 



|Gott|1972||Binney & Tremaine|1987] >. It follows that the majority 
of halo finders are based on the model and the free parameter that 
defines a given halo finder, e.g. the mean density of a halo or the 
linking length that defines friends-of-friends halos, is derived from 



the top-hat model. No such model exists that can provide a quanti- 
tative estimation of the value of one or more parameters that define 
the web. Hence, web classifier algorithms are bound to depend on 
some free parameters that cannot be determined from first princi- 
ples. The web that emerges from the geometric web classifiers de- 
pends strongly on the type or properties of objects used to trace the 
web, such as the mass of DM halos in simulations or the luminosity 
and spectral type of observed galaxies. The web that emerges from 
the dynamic web classification depends on the chosen smoothing, 
i.e resolution, and the value of the adopted threshold. No attempt is 
made here to optimize the values of these two parameters as the op- 
timization depends on the specific problem one wishes to address 
with the V-web. Different problems end up with different 'best' pa- 
rameters. We have chosen to use parameters that are best suited for 
the simulation studied here and the visual impression it leads to. 

The paper consists of a description of the algorithm that classi- 
fies the cosmic web (§2J, a brief summary of the N-body simulation 
used here (S|3} and a presentation of the results (*j4]). The paper con- 
cludes with a discussion of the V-web classification and its merits 
with respect to the analysis of structure and galaxy formation. 



2 WEB CLASSIFICATION: ALGORITHM 

Consider a DM-only N-body simulation. A description of the sim- 
ulation can be provided by the density and velocity fields evaluated 
on a finite grid, p(r) and v(r). The T-web is defined by the eigen- 
values of the Hessian of the gravitational potential which obeys the 
(rescaled) Poisson equation: 

V 2 (t> = A. (1) 

Here (f> is the gravitational potential rescaled by AnGp, p is the 
mean cosmological density, and the density (p) is given by A = 
pj p. The tidal tensor is defined as the Hessian of <f>, namely 



(2) 



dr a drp ' 

with a, ft — x,y and z. The three eigenvalues of the tidal tensor 
are denoted here by \J , where i = 1,2 and 3. (Throughout the 
paper it is assumed that the eigenvalues are arranged in decreasing 
order.) Note that 



A(r) = AF(r) + Aj(r) + A^(r). 



(3) 



The V-web is defined in terms of the shear tensor, which is rescaled 
here and is written as: 

^ 1 / dv a 



-'a/3 



(4) 



where Ho is the Hubble constant. The eigenvalues of T, a p are de- 
noted here as A^ (i = 1,2 and 3). 

The T-web and V-web are classified here by the number of 
the eigenvalues above a threshold, A^ and Ajh. We have not at- 
tempted to determine the values of these thresholds but rather to 
treat them as free parameters which define the corresponding webs. 
Our guiding principle is to adopt values which best reproduce the 
visual impression of the cosmic web. 



3 SIMULATION 

The web classification is applied here to a simulation performed un- 
der the framework of the Constrained Local Universe Simulations 
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(CLUES) projec"J whose aim is to perform cosmological simula- 
tions that reproduce the local large scale structure in the Universe 
as accurately as current observation admits. The simulation is per- 
formed within the ACDM cosmology, with cosmological param- 
eters consistent with a WMAP5 cosmology: a density parameter 
Q m = 0.28, a cosmological constant J1a = 0.72, a dimensionless 
Hubble parameter h = 0.70, a spectral index of primordial density 
perturbations n = 0.96 and a normalization as — 0.817. The con- 
strained nature of the simulation is ignored here and it is treated as 
a typical realization of the ACDMAVMAP5 cosmology. 

This is a pure DM simulation performed by the Gadget2 code 
within a computational box with side length Z/box = 64/i~ 1 Mpc 
using 1024 3 particles, corresponding to a particle mass of m p = 
1.839 x 10 7 /!. -1 Mq. The density and velocity fields are evaluated 
on a 256 3 grid, using Clouds in Cells (CIC) interpolation. Two CIC 
constructions are used here, the full computational box (BOX64) 
and an inner box of side length 8/i~ 1 Mpc (BOX8) placed at the 
center of the computational box. The CIC interpolation of the ve- 
locity field is performed by number weighting over the particles in 
adjacent cells. Differentiation of the velocity field is performed in 
Fourier space, using FFT The CIC procedure smoothes the field 
under consideration over two grid units, and therefore no aliasing 
is introduced by the FFT transform. 



4 RESULTS 

Figure [TJpresents the analysis of the LSS delineated by one of the 
principal planes of the simulation. The analysis is applied to the 
BOX64 CIC field with a grid resolution of 0.25/i _1 Mpc. Both the 
density and velocity fields have been Gaussian smoothed with a 
kernel of R s — 0.25/i _1 Mpc so as to suppress numerical artifacts 
on the grid scale, in particular the unphysical preferred directions 
induced by the Cartesian grid. The top panels show the (scaled) 
density field, A(r), in log (left) and linear (right) scales. The log 
scale map clearly distinguishes between the under- and over-dense 
regions. The middle panels show the V-web (left) and T-web (right). 
The threshold values used here are \Y h = 0.44 and A t T h = 0.7. The 
threshold values of the two webs have been chosen so as to obtain 
the maximal spatial resolution and to provide the best visual match 
to the appearance of the density field. 

The grey scaling of the web maps is: white (voids), light grey 
(sheets), dark grey (filaments) and black (knots). It should be noted 
that the maps shown here represent planar cuts, namely slices of 
one grid cell (0.25ft, _1 Mpc) thickness, through a three dimensional 
structures, hence the sheets appear as roughly ID filamentary ob- 
jects while the filaments appear as isolated spots, reminding the 
appearance of the knots. The lower-right panel shows the (nega- 
tive) of the divergence of the velocity fields, namely the sum of \f 
at each grid point. 

Inspection of the different plots of Figure[T]reveals a rich cos- 
mic web. A few immediate conclusions follow: 

(i) The V-web provides a much superior description of the 
web structure that is revealed by the mass distribution, compared 
with the T-web. This is best exemplified by the dense object at 
(X, Y, Z) « (-2, 12, 0)/i _1 Mpc. The V-web reveals the delicate 
structure around the object while the T-web yields an unresolved 
blob over there. In fact, the T-web shown in the figure has been 



1 http://www.clues-project.org/ 



constructed with an optimal choice of the smoothing and thresh- 
old parameters, so as to achieve the maximum spatial resolution. In 
particular, a reduction in the Gaussian smoothing does not lead to 
an improvement in the resolution of the T-web. 

(ii) The web permeates through the under-dense regions, where 
large extended almost empty regions are bisected by sheets. 

(iii) The sheets prevail mostly in the under-dense regions, com- 
pared with the filaments and knots which dominate the over-dense 
regions. 

(iv) Some of the sheets located in the under-dense regions ex- 
hibit a two caustics structure enclosing a void region. Namely, an 
inner planar part has a diverging velocity field (i.e., V ■ v ^ 0, is 
sandwiched between two caustic with a converging velocity field 
(V ■ v < 0). Such double sheets appear only in the under-dense 
regions. 

(v) The velocity divergence field follows very closely the den- 
sity field, with almost a complete match of the under-dense with the 
positive —V • v regions. Such a correspondence is the key for the 
success of the V-web in tracing the density field. Note that in the 
linear theory of gravitational instability under-dense points are al- 
ways associated with a divergent flow. It follows that the emergence 
of under-dense regions with converging flows is a pure non-linear 
dynamical effect. The present V-web algorithm classifies these ob- 
jects as sheets. 

Figures |5] and [5] present a zoom on the inner {Sh^Mpcf of 
the simulation, showing the density fields in linear and log scale 
and the V-web at a Gaussian smoothing of R s — 0.125/i _1 Mpc. 
The density and velocity fields have been CIC-ed on a 256 3 grid 
spanning the zoom region. The figures show parallel slices of 
Z = -3, -2, -1, 0, 1, 2 and 3/j _1 Mpc so as to depict the 3D 
structure. The color coding of the density (linear scale) field has 
been chosen so as to make the structure at the A(r) w 1 level to 
be more apparent. The log scale density maps provide a clear pre- 
sentation of the over-dense regions. The 3D structure that emerges 
from these parallel plots is of a filament that runs roughly parallel 
to the Z axis across the slices presented here. The width of the fila- 
ment depends somewhat on the value of the threshold used to define 
the web. The threshold also dictates whether the filament appears 
monolithic or if it breaks into smaller parallel filaments. That fil- 
ament is embedded within a curvy wall, i.e. sheet, of A(r) « 1 
that runs across the slices from Z = —3 to 3/i _1 Mpc. The wall 
leaves its trace on the parallel planes as a (partial) horseshoe lin- 
ear structure running from (X, Z) ~ (—0.5, — 4)ft _1 Mpc through 
~ (0, 0) towards « (4, 3)/i -1 Mpc. That wall is clearly classified 
by the V-web as a sheet, it is clearly visible in the (linear scale) den- 
sity maps. The wall borders a coherent under-dense region on the 
lower-right quarter of all the parallel planes. Now this under-dense 
region, that runs parallel to the Z axis, is not monolithic. It is fur- 
ther split by density ridges, forming some more smaller scale sheets 
that cut through the under-dense big blob creating smaller voids. 
These ridges are under-dense themselves, yet they are denser than 
the voids they surround. Some other coherent voids, sheets and fila- 
ments that run perpendicular to the Z = const, planes are apparent 
as well. 

The two mid-panels of Figure [T| which depict the V- and T- 
web, clearly show that the improved spatial resolution of the V- 
web compared with the T-web. The comparison is further extended 
in Figure [4] Again, lacking a quantitative measure of the quality of 
the construction of the cosmic web we resort here to a visual com- 
parison of the constructed cosmic web with the underlying density 
field. The plots show the XZ plane of the zoom box presented in 
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Figure 1. The normalized density field and the cosmic web, based on the CIC density field of the full computational box (BOX64) spanned on a 256 3 grid and 
Gaussian smoothed on the scale of 0.25fo~ 1 Mpc: a.The density field presented by log A (grey scale correspond to overdense and dashed contours to under- 
dense regions. The solid contour represents the mean density, (upper-left panel), b. The linear density field (represented by grey scale color map; upper-right 
panel), c. The velocity based cosmic web generated with = 0.44 made of voids (white), sheets (light grey), filaments (dark grey) and knots (black). Note 
that the map presents a planar cut through the cosmic web, hence sheets appear as long filaments and filaments as isolated compact regions (middle-left panel), 
d. The (gravitational) tidal cosmic web generated with \J h = 0.7 (same color coding for the cosmic web, middle-right panel), e. The multi-scale corrected 
V-web (see Section]?] bottom-left panel), f. The divergent of the velocity field: the grey scale shading corresponds to — V ■ v ^ (bottom-right panel). 



The cosmic web 5 




Figure 2. Zoom on the center of the computational box showing the density fields in linear and log scale and the V-web at a Gaussian smoothing of R B = 
0.125ft - 1 Mpc. Each row contains the V-web (left), log scale density (middle) and linear density field (right) maps. The upper, middle and bottom rows 
correspond to the SGZ = 3, 2 and l/i -1 Mpc planes, respectively. (The color coding and threshold values of Figurelllare followed here.) 



the series of the Z cuts of Figures [2] and [3] The sheet that almost 
coincides with the YZ principle plane is clearly manifested in the 
orthogonal XZ cut. The plots also show the filament that runs very 
close to the Z axis. The fine details, on sub-Megaparsec scale, of 
the density field are much better reproduced by the V-web than by 
the T-web. 



5 SOME STATISTICS AND A MULTI-SCALE 
APPROACH 

A close inspection of the the density and V-web maps (Figures 
|T] - |4j shows a very clear correlation of the V-web classification 
with the local density. There is clear correspondence between the 
web type with density, with voids, sheets, filaments and knots pro- 
gressing from the most under-dense to the densest environments, 
respectively. Yet, there is no one-to-one correspondence between 
the web type and density, as is manifested by Figure 5, whose left 
panel shows the probability distribution of grid cells for the differ- 
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Figure 3. Same as Figure[2], for SGZ = 0,-1, -2 and -3h _1 Mpc. 
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Figure 4. A comparison of the V-web and T-web of the XZ-plane of the zoom box (shown in Figures [2]and[3] using the same resolution and threshold values 
of Figure^ a.The density field presented by log A (upper-left panel), b. The linear density field (upper-right panel), c. The velocity based cosmic web (bottom 
-left panel), d. The (gravitational) tidal cosmic web (bottom-right panel). The color and contour coding of Figure[T|is followed here. 



ent web elements as a function of the fractional density, A. The 
plot shows that virtually at all density levels cells of any web type 
can be found. In particular, some very dense cells with A « 10 3 
are tagged as voids. A close inspection reveals that all such points 
are located at the very inner regions of massive halos, where the ve- 
locity field reflects the virial motions within virialized halos. This 
is a clear artifact of the over-resolution of the CIC grid which pen- 
etrates too deeply into massive halos. Only a very few cells are 
affected by this over-resolution and the global statistics is virtually 
unaffected. A similar effect is found also for some of the cells with 
density exceeding the virial mean density (A « 340 for the cos- 
mological parameters assumed here), which are classified as sheets 
and filaments, and are clearly located within the very inner regions 
of massive halos. The over-resoluton conjecture is easily verified 
by constructing a lower resolution V-web and confirming that the 
pathology disappears with the coarser grid. 

The following multi-scale procedure effectively removes the 



pathological misclassified cells. A low resolution V-web is con- 
structed, in addition to the finest desired web. All high resolu- 
tion cells defined as voids and have an over-desity exceeding unity, 
A ^ 1, are given the web attributes of the coarse web. Similarly, all 
high resolutions cells identified as sheets and filaments whose over 
density exceeds the virial density assume the web classification of 
the coarse grid. Both the high and law resolution webs are evalu- 
ated on the same 256 3 grid and the resolution is controlled by a 
Gaussian smoothing. The coarse grid is defined here by a Gaussian 
smoothing of R g — lh~ Mpc. The probability distribution of grid 
cells of the corrected web is shown by the right panel of Figure 5. 
The corrected V-web is presented by the bottom-left panel of Fig- 
ure^ From here on the V-web is assumed to be corrected by the 
multi-scale approach, unless it is otherwise explicitly stated. 

Table 1 presents the volume and mass filling factors of the 
V-web. The numbers in parentheses correspond to the V-web con- 
structed without the multi-scale correction. The numbers presented 
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web elements volume filling fraction mass filling fraction 



voids 


0.68 (0.69) 


0.13 


(0.15) 


sheets 


0.27 (0.26) 


0.36 


(0.37) 


filaments 


0.046 (0.046) 


0.34 


(0.37) 


knots 


0.0036 (0.0035) 


0.17 


(0.11) 



Table 1. The volume and mass filling factors of the various web elements. 
The filling factors obtained without the multi-scale correction are given in 
the parentheses. 



here clearly depend on the threshold level and resolution. Given 
that the V-web calculated here closely matches the visual appear- 
ance implied by the mass distribution, the following remarks are 
of interest. Voids and sheets occupy 95% of the simulated volume, 
and only less than 1% is occupied by knots. The volume and mass 
filling factors of the sheets and filaments are in close agreement 
with the findings of FR09, for a threshold of unity. The sheets are 
more dominant by volume and mass in the V-web compared with 
the T-web. 



6 DISCUSSION 

A novel algorithm for classifying the cosmic web has been intro- 
duced here. The algorithm extends the dynamic classification of 
Hahn et al (2007) and of FR09 by studying the eigenvalues of the 
(velocity) shear tensor at any given point rather than those of the 
(gravitational) tidal tensor. The web classification depends on the 
number of eigenvalues above a given threshold. Each point is clas- 
sified as either void (0 eigenvectors above the threshold), sheet (1), 
filament (2) or a knot (3). The value of the threshold is taken as a 
free parameter, determined so as to provide the best visual match to 
the observed LSS. As was pointed out in the Introduction there is no 
model which can make quantitative predictions of the cosmic web 
at the current epoch non-linear universe. This renders the various 
cosmic web finders to depend on one or more free parameters that 
cannot be derived from first principles. The V-web, defined by the 
threshold and by the adopted resolution, can serve as a platform for 
studying, among other things, the properties of galaxies and halos 
with respect to their environs. The choice of the parameters defin- 
ing the web needs to be optimized according to the problem at hand. 
The focus of the present paper is on the presentation of the method, 
choosing the threshold and resolution parameters so as to recover 
the visual impression that emerges from the LSS of the simulation. 
Detailed study and analysis of the cosmic web, its relation to the 
properties of DM halos and simulated galaxies is to be presented in 
a forthcoming series of papers. In particular, future studies will fo- 
cus on studying within the CLUES framework the Local Group in 
the context of the local cosmic web. A first step in that direction is 
the study of the orientation of the angular momentum of parent and 
sub-halos with respect to the cosmic web (Libeskind et al.|2012| l. 

Figures[T]|2]and|3]clearly show that the kinematic cosmic web 
(V-web) resolves much smaller structure than the gravitational cos- 
mic web (T-web). The two possible definitions of the web coincide 
on the large scales, where the linear regime prevails, but they depart 
with going to small scales where the non-linear dynamics manifests 
itself. Why does the velocity field trace the underlying cosmic web 
better? The following argument might provide a plausible explana- 
tion. It is known that the skeleton of the cosmic web has already 



been delineated by the initial conditions (Zeldovich 1970, Bond 
|et al.|1996) . It is also well know that the density field, and hence 
also the gravity field, evolves away from the linear regime faster 
than the velocity field (Kitaura et al. 201 1 , and references therein). 
This suggests that the velocity field retains a better memory of the 
initial conditions than the gravitational field, and therefore is ex- 
pected to trace the cosmic web better. The present calculations sup- 
port that assertion, yet the problem deserves a more in-depth study. 

The cosmic web permeates the large scale structure (Bond et al 
1996), covering all regimes from the very under-dense to the dens- 
est regions. Although the cosmic web correlates with the matter 
density, it cannot be solely described by the local density. For ex- 
ample, vast under-dense regions are bisected by planar structures 
(e.g. |Aragon-Calvo & Szalay||2012[ > that are clearly manifested 
by the matter distribution and are tagged as sheets by the V-web 
finder. The under-dense blobs, of A < 1.0, are surrounded by two- 
dimensional sheets with A w 1.0. Sheets contain quasi-linear par- 
tially collapsed one-dimensional filaments, which in turn contain 
highly non-linear compact knots. Figure 5 shows that the sheets are 
associated with the quasi-linear density field. The DM halo mass 
function has a strong dependence on the ambient density, and the 
mass function of halos residing in low density regions is heavily 
skewed towards low mass halos(Lemson & Kauffmann 1999, Hahn 
|et al.|20 07 l. Hence sheets are expected to be populated by low mass 
halos, and therefor to host predominantly faint galaxies. It follows 
that magnitude limited redshift surveys of galaxies which trace the 
distribution of the more luminous galaxies are expected to be bi- 
ased towards the luminous knots and filaments, leaving behind the 
dimmer sheets to be virtually undetected. As telescopes become 
more and more sensitive and begin to probe dimmer and dimmer 
galaxies, the internal structure of under dense cosmic web types 
(i.e. voids and sheets) will come into focus. 
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